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ABSTRACT 

The flux and line shape of the fine-structure transitions of Neil and Nem at 12.8 and 15.55 /.tm 
and of the forbidden transitions of O I A6300 are calculated for young stellar objects with a range of 
mass-loss rates and X-ray luminosities using the X-wind model of jets and the associated wide-angle 
winds. For moderate and high accretion rates, the calculated Ne II line luminosity is comparable to or 
much larger than produced in X-ray irradiated disk models. All of the line luminosities correlate well 
with the main parameter in the X-wind model, the mass-loss rate, and also with the assumed X-ray 
luminosity — and with one another. The line shapes of an approaching jet are broad and have strong 
blue-shifted peaks near the effective terminal velocity of the jet. They serve as a characteristic and 
testable aspect of jet production of the neon fine-structure lines and the Oi forbidden transitions. 
Subject headings: ISM: jets and outflows - ISM: kinematics and dynamics - stars: formation - stars: 
low-mass - X-rays: stars 
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1. INTRODUCTION 

More than 50 detections of the Neil 12.8 /im line have 
been made in young stellar objects (YSOs) by the Spitzer 
Space Telescope (Espaillat et al. 2007; Lahuis et al. 2007; 
Pascucci et al. 2007; Ratzka et al. 2007; Carr & Na- 
jita 2008) and by ground based telescopes (Herczeg et 
al. 2007; Najita et al. 2009; van Boekel et al. 2009; Pas- 
cucci & Sterzik 2009). The neon fine-structure lines were 
predicted by Glassgold et al. (2007; henceforth GNI07) 
to arise in disk atmospheres irradiated by X-rays which 
ionize neon and generate the warm and high-ionization 
conditions needed to excite the lines. However, the for- 
mation of even a single low-mass star occurs in a multi- 
component system consisting of a collapsing cloud core, 
a disk, an accretion funnel and outflows arising close to 
the star. Thus various authors have suggested that the 
Neil 12.8 pan. line might be generated elsewhere in the 
system and not just by the disk (e.g., Meijerink et al. 
2008, henceforth MGN08; Alexander 2008; van Boekel 
et al. 2009; Najita et al. 2009; Flaccomio et al. 2009; 
Guedel et al. 2010). 

Information on the origin of the Ne n emission can 
be obtained from the velocity resolved profiles of the 
12.8 /xm line. Najita et al. (2009) discussed four such 
measurements: TW Hya (Herczeg et al. 2007); T Tau 
(van Boekel et al. 2009); AA Tau, and GM Aur (Na- 
jita et al. 2009). TW Hya and GM Tau are transitional 
disks, the former seen almost face-on and the latter at 
an inclination of 54°, whereas AA Tau is a classical T 
Tauri star (TTS) with a large inclination angle (75°). 
Pascucci & Sterzik (2009) have recently reported four 
more line shape measurements: TW Hya, T Cha, Sz 73, 
and CS Cha; T Cha, and CS Cha are also transitional 
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disks. Najita et al. (2009) concluded that the line profiles 
for TW Hya, AA Tau, and GM Aur are consistent with 
the emission arising mainly from a gravitationally bound 
disk atmosphere. Pascucci & Sterzik (2009) measured 
small blue shifts in three transition disks which they 
interpreted as arising from low- velocity (~ lOkrns -1 ) 
photo-evaporative outflows, following Alexander (2008). 
Observations of T Tau by van Boekel et al. (2009) show 
that the Neil 12.8 /mi emission in this case arises from 
more than one location within a complicated and only 
partially revealed triplet stellar system. T Tau is one of 
several accreting sources with very strong Ne II emission 
(Guedel et al. 2008). Its total Ne II luminosity is 20 times 
larger than measured for many revealed TTSs observed 
by Spitzer. The line shape of Sz 73 shows mainly blue- 
shifted emission extending well beyond — 100 km s -1 , in- 
dicative of a strong stellar outflow. Neufeld et al. (2006) 
detected the Neil 12.8 /tm line in HH objects in the HH 
7-11 outflow. 

We show here that jets from YSOs, long known to emit 
the forbidden optical lines of heavy atomic ions, also gen- 
erate the mid-infrared transitions of Ne n and Ne III. Us- 
ing the X-wind model, we calculate the luminosities and 
shapes of these lines and also the forbidden O I A6300 
transitions for comparison. We obtain results for classi- 
cal TTSs (Class I and II YSOs) with varying mass-loss 
rates, and we show that active jets can dominate the 
observed flux of the Nell line. We also show that the 
line profiles have a distinctive shape that offers observa- 
tional opportunities to verify the origin of the emission 
in strong outflow sources and to test the predictions of 
the X-wind model of jets. 

2. METHODOLOGY 

The X-wind theory of low-mass star formation (Shu 
et al. 1994) describes how the interaction of the magne- 
tosphere and the accretion disk truncates the disk and 
forces its inner edge to co-rotate with the star. The in- 
teraction also drives an out-flowing wind and in-flowing 
accretion streams. In steady state, the inner edge of the 
disk coincides with the gravitational "X-point" of the 
system. Assuming that the flows emanate from this lo- 
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TABLE 1 

Parameters of the Reference Model 



Parameter 


Symbol 


Reference Case 


Stellar mass 


M* 


0.8 M 


Stellar radius 


fi, 


3.0 R Q 


Disk truncation radius 


fix 


4.8 R* 


Stellar rotation period 


2tt/Q x 


7.5 d 


Wind mass-loss rate 


Af„ 


1.6 X 10- 8 M Q yr- 1 


Disk accretion rate 


M A 


1.2 X 10- 7 M Q yr- 1 


Mean terminal wind velocity 


v w 


195 km s" 1 


Stellar luminosity 




2L @ 


X-Ray luminosity 


Lx 


3 X 10 31 ergs" 1 


Heating coefficient 




2 X 10~ 3 



cation allows a semi-analytic theory to be developed that 
yields a clear picture of outflows from low-mass YSOs. 
Shu et al. (1995) showed that the density asymptotes to 
cylindrical contours while the streamlines become radial. 
The physical conditions (temperature and ionization) of 
the X-wind were developed by Shang et al. (2002; hence- 
forth SGSL) using Shang's (1998) semi-analytic represen- 
tation of the flow. SGSL assumed the gas is atomic and 
heated by shocks and ionized by X-rays. They were able 
to give a good account of the optical observations of the 
forbidden lines and of the radio continuum observations 
of the Class I source, L1551 IRS 5 (Shang et al. 2004). 

In addition to X-rays, SGSL considered a variety of 
heating and ionization sources (listed in their Table 1). 
For example, in addition to X-rays, they treated H~ 
photo-detachment, Balmer continuum photoionization 
by stellar photons and by UV radiation from accretion- 
funnel hot spots, as well as collisional ionization. The 
X-rays can ionize a large part of the flow because radia- 
tive recombination occurs on a much longer time scale 
than the flow time scale (Bacciotti et al. 1995). Thus a 
significant level of ionization generated near the source 
of the flow is frozen into the wind out to larger distances. 
Many heating and cooling processes were considered, but 
SGSL found that the most important were: (1) adiabatic 
cooling of the expanding flow and (2) heating via the 
dissipation of internal velocity fluctuations and shocks. 
SGSL formulated the latter with a phenomenological for- 
mula for the mechanical heating rate per unit volume, 



where p and v are the local gas density and flow veloc- 
ity, s is the distance the fluid element has traveled along 
a streamline to the point of interest, and ah <C 1 is a 
phenomenological coefficient that characterizes the mag- 
nitude of the mechanical heating. Equation (TTJ) may be 
considered a prescription for shock heating in the jet. If 
velocity fluctuations of magnitude Sv generate the me- 
chanical heating, then a^ir ~ (Sv) 3 , and the values of 
ah needed to heat the flow, ~ 0.001, correspond to mod- 
erate velocity fluctuations, Sv/v ~ 0.1. Values of this 
order may be inferred from high-spatial resolution obser- 
vations of the optical forbidden lines emitted by revealed 
jets (e.g., Bacciotti et al. 2000; Woitas et al. 2002). 

We use the SGSL model to calculate the neon fine- 
structure line emission for the wind and jet of a low-mass 
YSO. We also compare the neon lines with Oi A6300, 
representative of optically forbidden lines. The results 



depend on the model parameters listed in Table 1. The 
most important are the wind mass-loss rate and X-ray 
luminosity In Table 1, both the wind mass-loss rate M w 
and the X-ray luminosity Lx are for a one-sided jet. The 
numerical values in the table are based on prior studies 
of jets (SGSL; Shang et al. 2004). We use the fiducial 
case of SGSL as our own reference model for a revealed T 
Tauri source that undergoes active accretion and drives 
a bright optical jet. 

SGSL conceived the X-ray emission to arise from a 
soft coronal source of dimension i?* enhanced by mag- 
netic reconnection regions, one located in the disk mid- 
plane and two others associated with magnetic field 
y-configurations above and below the plane (see the 
schematic drawing of Shu et al. 1997). They idealized 
this complex as a series of point sources, one at the ori- 
gin and the others displaced along the axis at heights 
±i? x with R x ~ 0.07 AU. The total X-ray luminosity 
of each source was composed equally of soft and hard 
X-rays, described respectively by thermal spectra with 
Tx = IkeV and Tx = 2keV. Because of the small di- 
mensions of all of the sources (~ Rx), they behave much 
like a single coronal X-ray source with an additional hard 
component. The soft X-rays are absorbed over relatively 
short distances for outflows with appreciable mass loss. 
According to the reference model, the more penetrating 
hard X-rays irradiate one hemisphere with a luminosity, 
L x ~ 1.5 x 10 31 ergs" 1 . 

The relatively large fiducial value of the X-ray luminos- 
ity in Table 1 was adopted by SGSL to treat active YSOs 
with bright jets and large accretion rates. Their concep- 
tion of the X-ray properties of YSOs was formed by the 
results from the early X-ray observatories, ASCA and 
ROSAT, as discussed by Shu et al. (1997) who related 
the properties of the fluctuating X-wind model and the 
emission of X-rays. Following Shu et al. (1997), SGSL 
based their X-ray luminosity on ROSAT observations of 
embedded sources (Neuhauser 1997), and took into ac- 
count the likely effect of absorbing material close to the 
YSO that reduces the fraction of X-rays that escape to 
be detected by an X-ray observatory. Potentially an even 
more important absorber than the wind is the system of 
accretion columns, which carry a larger mass flux (e.g., 
Alexander et al. 2004, 2005). Screening by circumstellar 
material may help explain the long-standing puzzle con- 
cerning the larger X-ray luminosities of weak-lined versus 
classical TTSs, as suggested earlier by Gahm (1980) and 
by Walter & Kuhi (1981). Gregory et al. (2007) have 
demonstrated this effect by calculating the propagation 
of X-rays through the accretion funnels of a magneto- 
spheric accretion model. 

As discussed above, the calculations reported here 
cover a wide range of X-ray luminosities and not just 
the fiducial level in Table 1: Lx spans 3 orders of magni- 
tude, from L x = 3 x 10 29 to L x = 3 x 10 32 ergs" 1 . This 
range may be compared with the X-ray luminosities in 
the Guedel et al. (2010) correlation study of the Nell 
line, where the measured X-ray luminosity varies over 
the range ~ 10 29 — 10 31 ergs -1 . As we have just pointed 
out, it is difficult to make a direct connection between our 
model parameter Lx and measured X-ray luminosities. 
However, our main current goal is to deduce the overall 
trend of line luminosities with X-ray luminosity and not 
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detailed agreement for any particular object. This point 
is reinforced by the fact that the numerical results are af- 
fected by other model parameters which are not varied. 
For example, the effect of X-ray ionization at the base of 
the wind is affected by the location of our model point 
sources. Thus, if the elevated source is moved up in z, a 
smaller Lx would be required to generate the same level 
of ionization. 

3. RESULTS 

Figure 1 shows the physical properties of the inner 
flow for the X-wind reference model in Table 1. The 
number density of atomic hydrogen is displayed with a 
color scale in Figure 1(a); the ordinate is the vertical 
distance z and the abscissa is the cylindrical distance w 
in units of AU. Beyond w = 5 AU, the density is ap- 
proximately collimated cylindrically. The velocity field 
(solid lines) reaches near-terminal velocity within a few 
AU from launch and tends to be radial at large distances. 
The occurrence of a cylindrically stratified density in an 
expansive flow is characteristic of the magnetically colli- 
mated jet in the X-wind model. 

The electron density is important for the production 
of the neon lines because the fine-structure transitions 
are excited by electron collisions. Figure 1(b) shows the 
electron density in color and the X-ray ionization rate C 
as solid black lines. Because the ionization tends to be 
frozen in, the electron density in the inner jet decreases 
more slowly with w than the ionization rate, which de- 
creases rapidly due to absorption and inverse-square di- 
lution of the X-rays. Both the electron density and the 
ionization rate manifest a conical shape similar to that 
often seen near the base of optical jets. Isotherms are also 
shown in Figure 1(b) as white lines. The magnitudes of 
the electron density and the temperature are suitable for 
exciting the neon lines. For example, the critical density 
of the Neil 12.8 /jm line is 4 — 5 x 10 5 cm -3 for temper- 
atures in the range 5,000 — 10, 000 K. In the very inner 
part of the jet, w < 2 AU, the upper level of the tran- 
sition is populated close to thermal equilibrium, but at 
larger distances the excitation is sub-thermal. For the 
reference case, most of the neon is ionized with roughly 
equal fractions in Ne + and Ne ++ ; Ne ++ gradually gives 
way to neutral neon at large distances along the inner- 
most core of the jet, while it continues to account for 
~ 50% of neon in the wide-angle wind. 

We have calculated the luminosity of the neon lines 
for a range of models in which the main parameters, the 
X-ray luminosity Lx and the mass-loss rate M w both 
for a one-sided jet are varied away from the reference 
values in Table 1. In some cases, their ratio is kept fixed 
at Lx/M w = 3 x 10 13 ergg _1 , while the wind mass-loss 
rate changes from 5 x 10~ 7 to 5 x 10~ 10 Af©yr _1 . This 
range extends from a somewhat embedded Class I source 
down to the level of a revealed classical TTS. Line fluxes 
are calculated for the standard sequence of parameters 
used by Shang et al. (2004) . Then the mass- loss rates are 
increased and decreased by factors of 3 keeping Lx/M^ 
fixed. Finally, the constraint of fixed Lx/M w is removed 
for each of these standard cases, and Lx is varied by 1 
order of magnitude above and below the standard series. 

The ionization balance of the neon species is calculated 
by the method given in Section 2 of GNI07. Excitation 



by neutral hydrogen has been ignored because the large 
electron fractions make electron collisions dominant. The 
total abundance of neon is set at 10~ 4 relative to hy- 
drogen and that of oxygen at 6.6 x 10 -4 , following the 
discussion in GNI07. The populations of the lowest five 
levels of Ne ++ and O i were solved exactly as were the 
two levels for Ne + . The lines were assumed to be opti- 
cally thin, a good approximation for the region shown in 
Figure 1. 

Figure 2 plots the luminosity of the Neil 12.8 /im and 
Nelll 15.55 /im lines and the forbidden transitions of O I 
A6300 against three variables: (a) the mass-loss rate M w , 
(b) the X-ray luminosity Lx, and (c) the product LxM w 
- introduced by Guedel et al. (2008) in a search for 
empirical correlations of the Ne II luminosity with other 
properties of YSOs. Values from both groups of calcu- 
lations (fixed and variable Lx/M w ) are plotted. The 
fluctuations in Figures 2(a) and (b) show the effects of 
the variable not plotted. The quantity M w is intrinsic 
to X-wind theory, whereas Lx is an external parameter 
that was added in the development of the physical prop- 
erties of the X-wind by SGSL. There is no self-consistent 
theory available yet that treats both the magnetospheric 
wind and accretion generation together with X-ray emis- 
sion. The heating coefficient an in Equation [T] is another 
external parameter introduced by SGSL, but it is less 
effective in changing the physical properties of the out- 
flow, and for purposes of simplicity it is not varied very 
far from the value in Table 1 for the reference model. The 
large open triangles in Figure 2 are for the reference case. 
The model calculations all fit power laws with relatively 
small scatter, i.e., the 3u uncertainties in the least-square 
slopes are all significantly smaller than the slopes them- 
selves. The parameter ah has been kept fixed at 0.002 in 
Figure 2; varying ah by a factor of a few would naturally 
introduce some scatter without changing the slope. Un- 
like the case of disk emission (MGN08), the correlation 
with X-rays in Figure 2(b) is reasonably well founded 
because the effects of the other main parameter M w are 
seen in Figure 2(b) to be relatively modest|3 

Somewhat surprisingly, the fits with the empirical pa- 
rameter LxM w in Figure 2 are superior to the other two. 
The ratio of the ler uncertainty in slope to the slope itself 
is ~ 0.05. A possible interpretation of this correlation is 
that it reflects the primary role of the electron density in 
producing the neon ions and in exciting their line emis- 
sion. This in turn follows from the idea that Lx/M w is 
the global average of the local X-ray radiation parame- 
ter, Q/n. This quantity is roughly proportional to the 
square of the electron fraction, so that the square root 
of the empirical parameter, LxM w , is a measure of the 
electron density. 

The slope of the correlations for O I in Figure 2 is larger 
than those of the Ne n and Ne in lines, and the Ne II slope 
may be slightly larger than the one for Ne III. These dif- 
ferences can be traced to the ionization fractions of these 

4 Figure 2(b) is similar to Figure 16 of MGN08, where the lumi- 
nosity of several diagnostic lines, including the Nell 12.8 fim line, 
are plotted versus Lx f° r the same disk structure, that of a generic 
TTS disk. A correlation diagram for a real sample of observed TTS 
disks may not follow this power law because it will likely include a 
variety of disks at different ages and density structures. 
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Fig. 1. — (a) Background density of atomic hydrogen (color scale) and velocity contours (solid black lines) for the X-wind reference case, 
(b) Electron density (color scale), ionization rate £ (solid black lines) and temperature (white lines). Volumetric densities are given in cgs 
units and velocities in kms -1 . 



ions. As pointed out by MGN08, direct ionization of neu- 
tral oxygen by X-rays is not competitive compared with 
charge exchange with H + . As a result, the neutral oxygen 
fraction in the present calculations is above 90% except 
for very small values of Lx and M w . On the other hand, 
the Ne + and Ne ++ fractions are sensitive to the local 
X-ray ionization parameter, C/ n - Thus for the reference 
case, the Ne ++ /Ne + ratio approaches unity in the inner 
jet, but then falls to ~ 0.5 throughout the rest of inner 
region of the outflow (zu < 10 AU) due to attenuation 
and inverse-square dilution of the X-ray flux. For larger 
values of M w , attenuation takes an increasing toll and re- 
duces the Ne ++ /Ne + ratio, although it is always larger in 
the inner jet than in the outflow at large. This explains 
the tendency of the luminosity ratio L(Ne III) /L (Nell) 
in Figures 2(a) and (c) to decrease with increasing M w . 
This tendency is at most moderately significant because 
of the large fluctuations in the theoretical predictions of 
the luminosities of the Neil and Nem lines in Figure 2. 

The fact that the luminosities in Figure 2 correlate well 
with the variables M w and Lx means that they correlate 
with one another. This is illustrated in Figure 3, where 
the Neil 12.8 /jm luminosity is plotted against the Oi 
A6300 luminosity on a log-log scale. The slope of the 
correlation line is 0.83 ± 0.048. The O I A6300 lines have 
long been used as diagnostics of YSO jets. The predicted 
correlation with Neil suggests that the latter line may 
also be a good diagnostic for these jets. Because our 
calculations are based on X-wind theory supplemented 
by chemical-physical considerations, including X-ray ion- 
ization of the outflow, Figure 3 offers a potentially use- 
ful observational test of the excitation theory in SGSL 
without the uncertainties of external parameters such as 
Lx- A preliminary inspection of the available data for 
jet sources in Guedel et al. (2010) suggests a good corre- 
lation with our theory. Further detailed comparisons of 
the Neil 12.8 /mi and the Oi A6300 luminosities should 
be helpful in understanding the origin of these emission 
lines. 

Figure 4(a) shows normalized profiles of the Neil 12.8 
/xm, Neili 15.55 /jm and the Oi A6300 lines for the refer- 
ence case of Table 1 at an inclination angle of 45°. The 



calculations are for one side of the outflow and the results 
have been binned to a spectral resolution of 1 kms -1 . 
The predicted emission extends over a broad range of 
velocities from ~ —200 to ~ +100 kms -1 at this inclina- 
tion, The dominant feature is the large blue-shifted peak 
close to the projected terminal velocity of the jet. Fig- 
ure 4(b) shows the contributions to the Neil line shape 
as a function of the cylindrical coordinate zu. As zu in- 
creases, the height of the peak grows from contributions 
of regions with ever-increasing radius and outflow veloc- 
ities close to the terminal value. The Nell line has a 
broader low-velocity wing than the other lines, indicat- 
ing that Neil fraction is more sensitive to the physical 
conditions close to the source and the acceleration zone. 
On a larger scale, the Ne II line is generated at relatively 
small distances, zu < 5 AU, with the broad wing from 
~ —50 to ~ +50 kms" 1 produced within ~ 1 — 2 AU of 
the axis. This is close to the source of the wind; Fig- 
ure 1(a) shows that the wind velocity is already quite 
large by z = 0.5 AU. 

The line profiles vary with the inclination angle, which 
is important for making comparisons with observations. 
Figures 5(a) and (b) show line shapes at 30° and 60° for 
one side of the jet; they are similar to the 45° case shown 
in Figure 4(a). Significant changes occur for inclinations 
near 90° and 0°. Near these extremes, the emission peak 
occurs close to the real terminal velocity for the face-on 
case (0°), while it approaches the stellar velocity for the 
edge-on case (90°). The X-wind jet flow has the form of a 
cylindrically symmetric rounded cone or bowl. The peak 
velocity in the emission is determined by the projection 
of the terminal velocity. The line shape varies gradu- 
ally with inclination angle as the projected velocity peak 
moves toward larger or smaller projected values. The 
low velocity wing remains broad and extends to negative 
velocities. In theory, the jet is bipolar, and the full pro- 
files are in principle symmetric with respect to the stellar 
velocity. The receding jet produces a red-shifted emis- 
sion peak at the same magnitude of projected velocity as 
the blue-shifted peak, and its low- velocity wing extends 
to the blue. Depending on the actual viewing angle and 
orientation of the system, the line profile of a jet with 
two visible lobes will be a mixture of shapes from each 
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Fig. 2. — (a) Line luminosities Nell 12.8 /im (solid line), Neoi 
15.55 /an (dotted line), and the Ol A6300 (dashed line) from a 
one-sided jet, plotted versus M„ in units of Mq yr — . The vertical 
axis labels line luminosities in erg s — 1 . The open triangles mark 
the reference case of Table 1. Af w is the mass- loss rate from one 
hemisphere. The slopes of the power-law fits are: 1.351 ± 0.132 
(Ol); 1.215±0.084 (Nell); 0.851±0.110 (Nem). (b) Same as (a), 
except the line luminosities are plotted versus the X-Ray luminosity 
Lx of a one-sided jet in units of erg s — 1 . The slopes of the power- 
law fits are: 1.328 ±0.117 (Ol); 1.001 ±0.152 (Nell); 0.813 ±0.068 
(Ncni). (c) Same as (a), except the line luminosities are plotted 
versus LxAf w in units of erg s^ 1 X Mqjt -1 . The slopes of the 
power-law fits are: 0.709±0.022 (Ol); 0.606±0.026 (Nell); 0.468± 
0.022 (Ncin). 
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Fig. 3. — Luminosity of the Nell 12.8 (im line vs. the Ol A6300 
luminosity in a log-log plot; the units are erg s — 1 . The correlation 
line has a slope close to 0.8. 

side. If for some reason only the receding side is visible, 
the profile will be peaked at the high red-shifted veloc- 
ity. Several authors have suggested that sources with 
observed high Neil luminosity are associated with jets 
(e.g., Guedel et al. 2010; van Boekel et al. 2009; Flac- 
comio et al. 2009). Our results support this view, and 
our specific predictions can be checked by the detection 
of strong blue-shifted line emission, as seen in Figures 4 
and 5. 

4. DISCUSSION AND CONCLUSIONS 

We have calculated the emission of the Neil 12.8 /im, 
Nem 15.55 fim and Ol A6300 lines for the X-wind model 
for the outflows from YSOs. The oxygen lines are well- 
established diagnostics of jets, whereas the Nell 12.8 /im 
line has only recently become a probe of the circumstellar 
gas of young stars, as described in the Introduction. We 
have calculated the dependence of the line emission on 
the two main parameters of the model, the wind mass- 
loss rate M w and the X-ray luminosity Lx- The mass- loss 
rate is an intrinsic parameter of X-wind theory, whereas 
the X-ray luminosity was added by SGSL in modeling the 
physical properties of the wind. Figure 2(a) expresses 
our first basic result. It shows that the line luminosi- 
ties correlate well with the X-wind parameter M w . Fig- 
ure 2(a) also shows that, for mass-loss rates in excess of 
those characteristic of Class II TTSs, the luminosity of 
the Nell 12.8 /im is well above the observed range for 
this line given by Guedel et al. (2010). It is interest- 
ing that the combination parameter, LxM w , gives the 
best correlation of the three shown in Figure 2. Figure 3 
gives another example of a good correlation, in this case 
between the O I forbidden transitions and the Ne II fine- 
structure line. It suggests that the Neil 12.8 /xm line is 
an excellent probe of the jets from YSOs, as are of course 
the well-established O I A6300 line. 

The theoretical line profiles displayed in Figures 4 and 
5 have a distinctive shape. For an approaching one-sided 
jet, they have a strong blue-shifted peak near the wind 
terminal velocity that is associated mainly with the inner 
part of the jet within w < 5 AU. There is also a broad 
shoulder centered around the stellar velocity that arises 
near the source of the jet and the acceleration region of 
the outflow. These features reflect some of the unique 
properties of the X-wind model of outflows from YSOs. 
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Fig. 4. — (a) Line profiles for Nell 12.8 /im (blue), Nelll 15.55 /im line (red) and the Ol A6300 (black) for the X-wind jet, calculated for 
the reference case of Table 1. The inclination angle is 45°, and the line shapes are normalized to their peak values, (b) Buildup of the Nell 
line shape as a function of cylindrical radius vj for the case of 45° inclination angle. 

irradiated disk, either from the disk proper (MGN08) or 
from a photo-evaporated wind (Owen et al. 2010; Er- 
colano & Owen 2010). 

Line profile measurements have already been carried 
out for seven low-mass YSOs, as discussed in the Intro- 
duction. Four of the seven are transition objects which 
have evolved disks with substantial gaps in the inner 
dust distribution (TW Hya, GM Aur, CS Cha, and T 
Cha), typically on the AU scale. These four objects all 
show clear signatures of accretion onto the stellar sur- 
face, which implies that there is a significant amount 
of gas present within the inner radius of the dust gap. 
These are not the most appropriate systems for apply- 
ing the X-wind model, originally conceived to describe 
YSOs with strong outflows. Two of the seven YSOs (T 
Tau and Sz 73) are outflow sources. The simple X-wind 
model used in this paper is inadequate to deal with the 
T Tau triple stellar system, which has multiple sources 
of X-rays and Ne II line emission. Of all the line shapes 
reported by Pascucci & Sterzik (2009), the broad, blue- 
shifted line of Sz 73 is closest to the profiles we show in 
Figures 4 and 5. However, little is known about this 
YSO, and it would be premature to conclude that it 
agrees (or disagrees) with our model until more infor- 
mation is obtained, e.g, its X-ray properties determined, 
and detailed modeling is carried out. This leaves AA 
Tau, whose line shape has been measured by Najita et 
al. (2009) to be broad (FWHM - 70 km s" 1 ) and with 
a red-shifted peak at -rT5kms _1 . There is also emis- 
sion corresponding to the blue-shifted velocity, but the 
spectrum is noisy and difficult to characterize. With its 
large inclination angle, it is possible that higher signal- 
to-noise observations may show the line shape expected 
on the basis of our X-wind calculations, especially since 
AA Tau manifests some evidence for outflow (Bouvier et 
al. 2003, 2007). On the other hand, the flux of its Neil 
12.8 fxm line is at the level obtained for X-ray irradiated 
disk models (MGN08). Thus any serious model calcu- 
lation would have to include all the sources of emission, 
the accretion funnel, the disk and the wind. We would 
like to recommend that searches be made for Nell 12.8 
/im line emission from strong jets already detected in for- 
bidden transitions from O I and other heavy ions using 
high-resolution MIR spectrometers on large telescopes. 

Hollenbach & Gorti (2009) have provided an alterna- 
tive view of the Neil 12.8 /im line emission for sources 




-200 



-100 
V [km s - 

(b) 



1 00 200 



Fig. 5. — (a) Line profiles for Nell 12.8 /im (blue), Nelll 15.55 /im 
line (red) and the O I A6300 (black) for the approaching X-wind jet, 
calculated for the reference case of Table 1. The inclination angle 
is 30°, and the line shapes are normalized to their peak values, (b) 
Line profiles for Nell 12.8 /im (blue), Nelll 15.55 /im line (red) 
and the Ol A6300 (black) for the X-wind jet, calculated for the 
reference case of Table 1. The inclination angle is 60°, and the line 
shapes are normalized to their peak values. 

Observations to test these predictions for high luminos- 
ity sources are feasible for both the Neil 12.8 /im and the 
O I A6300 lines. High-spectral resolution measurements 
would be of great interest for understanding the source 
of the Ne II and O I line emission, the role of X-rays in 
determining the physical properties of outflows, and the 
dynamics of jet formation. The blue-shifted peak of the 
Neil line is probably the most important signature of 
jet emission, because the low-velocity shoulder may be 
contaminated in some cases by emission from the X-ray 
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with high accretion rates. In their theory, the main 
source of ionization of neon and the excitation of the 12.8 
pm line is high-velocity (~ 100 kms -1 ) shocks that affect 
the entire jet. While high-velocity shocks may generate 
large enough temperatures to ionize neon in localized re- 
gions, there is little observational evidence for the strong 
and continuous shocking of jets. The bulk of the mag- 
netized flow is at most mildly shocked, despite the high 
terminal speeds reached close to the source. It is these 
mild shocks that are responsible for heating the outflow 
in the present model. Hollenbach & Gorti (2009) do not 
calculate line shapes for their high- velocity shock model. 
Although they are likely to be broad (<~ 100 km s -1 ), the 
line shapes can be expected to be different from the ones 
shown here because the broadening arises downstream 
along the jet where the strong shocks occur. In con- 
trast, most of the Ne n and Ne in emission in our work is 
generated within or close to the acceleration region that 
surrounds the source, so that the line shape takes its final 
form close to the source, as shown in Figure 4. 

Guedel ct al. (2010) have made extensive searches for 
observational correlations of the Nell 12.8 pm line lu- 
minosity with X-ray luminosity and other properties of 
YSOs, as have Flaccomio et al. (2009). They do not find 
any particular trend for low-accretion systems, where 
the measured Neil luminosities span a range of a fac- 
tor of 10 or more about a median of approximately 
L(Neii) ~ 5 x 10- 6 L Q (or log i (Nell) = 28.25 in cgs 
units). This value is close to our reference theoretical 
value, the large open triangle in Figure 2. It is also sim- 
ilar to the value predicted by MGN08 for an X-ray ir- 
radiated disk of a typical TTS, although the Ne n emis- 
sion from disks has not been calculated for any other 
case. Calculations for photo-evaporated winds give sim- 
ilar values (Alexander 2008; Ercolano & Owen 2010). 
One explanation for the large number of measurements 
for both optically thick and transition disks with this 
level of Ne II luminosity is that other variables are im- 
portant, and not just the X-ray luminosity. Schisano 
et al. (2010) illustrate this possibility by showing that 
variations in disk flaring and X-ray spectrum give rise 
to a scatter in the predicted Neil luminosity similar to 
what is observed. According to our calculations, outflows 
may also contribute in some cases, but they are not the 
dominant source of emission for low or even moderate ac- 
cretion rates and X-ray luminosities. However, as shown 
in Figure 2, they may dominate for high mass- loss rates 
and X-ray luminosities. Thus our calculations support 
the "bi-modal" picture of Guedel et al. (2010) and van 
Boekel et al. (2009), especially their suggestion that the 
Nell luminosity of high- accretion sources arises in out- 
flows. 

On the basis of a regression analysis, Guedel et 
al. (2010) also concluded that the Neil luminosities of 
jet sources correlate, not only with X-ray luminosity, but 
with accretion rate and mass-loss rate. The significance 
of these correlations is diminished by the relatively small 
sample of sources for which the parameters have been 
measured and by scatter in the data. This is especially 
true in the case of the mass-loss rate, where only two 
sources (T Tau N and DG Tau) have exceptionally large 
Nell luminosities. Without relying on these two "out- 



liers" the correlation of the Ne n luminosity with mass- 
loss rate is not yet definitively established. Such a corre- 
lation would be directly relevant to the present calcula- 
tions because, in X-wind theory, it is the mass-loss rate 
that directly determines the properties of the outflow, 
rather than the accretion rate, although the two are re- 
lated. Mass-loss rates of YSO jets are often determined 
from bright optical emission lines such as O I A6300 (e.g., 
Hartigan et al. 1994, 1995), but the results for YSOs of 
the same mass can range over several dex (White & Hil- 
lenbrand 2004). In light of such uncertainties, it would 
be better to compare the luminosity of the O I A6300 with 
that of the Nell line, as we do in Figure 3. Guedel et 
al. (2010) test this correlation with data for sources they 
identify as optically thick disks without jets, transition 
disks, and jet sources. A preliminary inspection of the 
available data for the seven jet sources in their Figure 4 
with L(Nell) > 8 x 10~ 6 L© suggests a good correlation 
with our theory. Further observations of the Nell and 
O I emission from strong jets including line shape mea- 
surements would be extremely useful in this connection. 

The correlation predictions in Figure 2 suggest an- 
other possible test of our calculations using the Nelii 
15.55 pm line, although only a few detections of this 
line have been reported so far. In Figure 2(a) the pre- 
dicted but uncertain ratio of the Neil 12.8 /im to the 
Nelii 15.55 /xm line flux is in the range ~ 1 — 10; the 
Neni flux is ~ 10~ 14 — 10 _13 ergcm -2 s _1 . Lahuis et 
al. (2007) reported a tentative detection of the Ne in line 
with Spitzer for Sz 102, a poorly studied source with 
a prominent jet, known to emit primarily soft X-rays 
(Guedel et al. 2010). Lahuis et al. (2007) give the neon 
line fluxes as L(Neii) = 3.6 x 10 -14 ergcm -2 s _1 and 
L(Neiii) = 2.3 x 10 _15 ergcm _2 s _1 , or a Nelll/Nell 
luminosity ratio of ~ 1/16. This small value is consis- 
tent with a disk origin for the neon lines, but the mea- 
sured fluxes are smaller than predicted by MGN08 for 
the generic TTS disk and the measured X-ray luminos- 
ity. Flaccomio et al. (2009) detected both neon lines from 
the F7 Class III YSO WL5/GY246 in the p Oph clus- 
ter. Although the luminosities are in rough accord with 
a disk origin, the Nem/Neii luminosity ratio of ~ 1/4 
is somewhat higher than given by MGN08 and more like 
the values obtained here for jets. However, we would not 
expect to find a strong jet in a Class III YSO, so the 
disk model is favored in this case. A careful determina- 
tion of upper limits if not detections in Spitzer spectra of 
high luminosity sources of the Nelii 15.55 pm line could 
provide further checks of our theory. Since Spitzer opera- 
tions have entered the warm phase, further direct detec- 
tions of this line from space will not be available for some 
time. However the optical lines of Nelii A3869/3967 are 
accessible from the ground, and they have been observed 
from Hn regions and other highly-ionized sources. Pre- 
liminary calculations indicate that these lines have about 
the same luminosity as the Neni 15.55 pm line. We sug- 
gest that these lines be used to further probe the ioniza- 
tion state and emission characteristics of jets and disks 
around low-mass YSOs. 

In conclusion, we have calculated the emission of the 
fine-structure lines of neon and of the forbidden O I tran- 

5 T Tau N is an outlier in that it is part of a compact triple 
system; DG Tau has an unusually small observed X-ray luminosity. 



sitions near 6300 A, according to the X-wind model of 
jets. These lines trace similar aspects of the outflow. 
They correlate well with the X-wind model parameters 
and with the combination variable LxM w , as well as with 
one another. We have suggested several tests of the cal- 
culations, including the measurement of the unique blue- 
shifted peak in the line profile near the wind terminal 
velocity. We also suggest that jets may account for the 
luminous mode in the bimodal picture of Ne n correla- 
tions proposed by Guedel et al. (2010) and van Boekel 
et al. (2009). 
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